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Abstract 

The energy evolution of a quantum chaotic system under the perturbation that harmonically 
depends on time is studied for the case of large perturbation, in which the rate of transition 
calculated from the Fermi golden rule exceeds the frequency of perturbation. It is shown that the 
energy evolution retains its diffusive character, with the diffusion coefficient that is asymptotically 
proportional to the magnitude of perturbation and to the square root of the density of states. The 
results are supported by numerical calculation. They imply the absence of the quantum-classical 
correspondence for the energy diffusion and the energy absorption in the classical limit 7l — ^ 0. 
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I. INTRODUCTION 



The problem of susceptibility of chaotic systems to perturbations has attracted much 
attention in the last decade [1 - 9]. This problem is fundamental, since it includes the 
determination of the response of a material system to an imposed external electromagnetic 
field, the setup that is typical for many experiments. Due to the sensitivity of classical phase 
trajectories or quantum energy spectra and stationary wavefunctions of chaotic systems to 
small changes of their parameters, the problem is challengingly difficult. A consistent and 
noncontroversial picture covering (albeit qualitatively) all the essential cases of the response 
hasn't been yet drawn at present. From the point of view of general theory, the problem 



is related to the applicability of the concept of quant um-classica 
systems, that is a long-standing question in its own right 
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correspondence to chaotic 



We shall study a one-particle system with the Hamiltonian of the form H = H — 
Fxcosuiot, where Hq (p, f) is the Hamiltonian of the unperturbed system; p and f are the 
operators of Cartesian components of the momentum and of the position of the particle. 
The classical system with the Hamiltonian function H (p, r) will be assumed to be strongly 
chaotic, that is, nearly ergodic on the energy surfaces in a wide range of the energy values, 
system with d > 2 degrees of freedom. In the perturbation operator V (t) = —Fx coswo^ the 
active variable x is one of the Cartesian coordinates of the particle, coupled to the external 
homogeneous force field. The amplitude F in the following will be referred to as field. In 
the following we shall deal with the quasiclassical case, when the Planck constant is small 
in comparison of the action scale of the system H . 

Under the influence of the perturbation the energy value E(t) = H (t) varies in a quasir- 
andom way. These variations frequently can be described as a process of the energy dif- 



fusion 
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131 ] . when for the ensemble with the microcanonical initial energy distribution 
H (0) = E the dispersion of the energy increases linearly with time, (AE 2 (t)) = 2Dt, where 
D (E, F, uq) is the energy diffusion coefficient. 

If the external field F is sufficiently small in comparison with the appropriately averaged 
values of the forces acting on a particle in the unperturbed system, then in the classical 
model the energy diffusion coefficient D can be expressed through the characteristics of the 
unperturbed chaotic motion of the active coordinate, namely 

D = ^u 2 F 2 S x (E,u ), (1) 

where S X (E, uq) is the power spectrum of the active coordinate (the Fourier transform of its 
autocorrelation function) for the motion on the surface with the constant energy value E [9] . 



The same expression (1) in the case of weak perturbation can be obtained in the classical 
limit from the quantum model. The evolution of the quantum system can be treated as 
a sequence of one-photon transitions between stationary states of the unperturbed system 
|n) — ► | A;), accompanied with the absorption or emission of the quanta Hluq. For small h the 
energy spectrum of Hq is quasicontinuous, thence the rates of transition are given by the 
Fermi golden rule (FGR) 

W F = ^-F 2 \x nk \ 2 p(E k ), (2) 

where x nk is the matrix element of the active coordinate, and p{E k ) is the density of states 
near the final state of the transition. Although the matrix elements x nk in quantum chaotic 
systems fluctuate wildly with the variation of k [10, 11], the averaged squared quantity 
\x n k\ 2 in the limit h — > is smooth; it is proportional to the power spectrum S x (E,u>o) of 
the coordinate 



14. 15], 



S x (E, LO ) 



lXnhl ~ hp (E) ' { ) 

From Eqs. (2) and (3) we have for the transition rate 

W F = ^F 2 S x (E,u). (4) 

Then for the energy dispersion for small t we have (AE 2 ) = 2 (Hooq) 2 Wpt, that returns us 
to the Eq.(l) for the energy diffusion coefficient. It can be shown that the same expression 
for D holds also for large t The energy absorption in chaotic systems comes as an 
epiphenomenon of the energy diffusion [4]. With the account of the dependence on the 
energy of the power spectrum S x (E, uj) and the density of states p (EYthe diffusion becomes 
biased, and the energy absorption rate Q is given by the formula 3, ^ 

Q-~ipD). (5) 

Although for weak fields D does not depend on the Planck constant h, the condition of 
the applicability of Eq. (2) does. The FGR is, after all, only a formula of the first order 
perturbation theory. It is based on the assumption that the transition process has a resonant 
character - that the width A of the energy distribution of states populated from the original 
one, given by the Weisskopf - Wigner formula A = hW |l6j], is small in comparison with 
the quanta energy Huq. From Eq. (4) it is evident that in the classical limit h — > this 
condition will be violated. In the following we shall use the border value of the field Fb, 
defined by the condition Wp (F b ) = u , and refer to the domain F > F b as the range of the 
strong field. 



By analogy with other models, for strong fields one can expect a slow-down of the growth 
of the energy diffusion coefficient D and of the energy absorption rate Q. For example, for a 
two-level system with relaxation the quadratic dependence the absorption rate Q oc F 2 for 
small field turns into a field-independent value Qo for strong one. The border is determined 
by the condition Q 2 /TiT2 ~ 1, where Q is the Rabi frequency and I^r^ are longitudinal 



and transversal relaxation rates correspondingly 17]. The rate of transitions from discrete 
to continuous energy spectrum (that are basically covariant with the energy absorption rate 
Q), studied in the context of the theory of photoionization, for sufficiently strong fields can 
even decrease with the increase of F - the effect that is known as atom stabilization by the 
strong field [3| . 

For our model we can describe the slow-down by a simple estimate. Assuming the expo- 
nential decay of the initially populated state and using the Weisskopf - Wigner formula for 
the energy distribution of the final states [16], we can rewrite the expression for the rate of 
the transitions (with the absorption of quanta) in the form 

1 W 



w 71 rV [bJ] 



2 J h 



7T (uj - LOof + W 2 



p (u) du, (6) 



where V(u) is a matrix element of perturbation taken as a function of the transition fre- 
quency and uj = (Ek — E n )/h. In the weak field the expression in square brackets turns into 
the 5 - function, and Eq. (6) returns the usual FGR. However, if we dare to use Eq. (6) as 
an equation that is valid for any magnitude of the perturbation, then in the case of strong 
fields we will have 

W 2 = \j^ P {u)du. (7) 

The rate of transitions W, and, consequently, the energy diffusion coefficient D and the 
energy absorption rate Q, become proportional to the magnitude of perturbation F and to 
the square root of the density of states. This derivation of Eq. (7) is too risky - but in the 
following sections this conclusion will be put on a more solid footing. 

The slow-down of the energy diffusion in quantum chaotic systems in strong harmonic 
fields has been first demonstrated by Cohen and Kottos [5]. However, their analytical 
estimates and data of numerical experiments are in quantitative disagreement with the 
results of the present paper. 

Lastly, it must be noted that the strong field regime is easily attainable to experiments. 
For example, for the excitation of multiatomic molecules with the infrared laser radiation 
the border field corresponds to the intensity value / ~ 10 9 W cm" 2 , that has been reached 
in experiments long time ago jlo| . 



II. THEORY 



For the system with the Hamiltonian H = Hq + V cos uot we take the wave function in 
the form of the expansion in the basis of stationary {(p m } states of Hq, 

^(t) = Y.a m Vm(v)e- l ^ t . (8) 



i-JT = J2 n km cos uj te iuJkmt a m , (9) 



Then for the amplitudes {a m (t)} we obtain the system of equations 

da i 

dt k 

where the quantities Vt^n = Ti~ l Fxk n are the Rabi frequencies of transitions. We shall use 
the initial conditions a m (0) = 8rn n \ at the initial moment only one of stationary states, 
tp n , is populated. Following [l^, U^, we shall assume that x n k are independent random 
Gaussian variables with zero mean value and the dispersion, given by Eq. (3). The system 
of equations (9) will be treated as a member of the corresponding statistical ensemble. 

We are going to concentrate on the process of the energy diffusion. Then in the zeroth 
approximation we can restrict ourselves by consideration of the evolution of probability 
density in a narrow energy range around the initial state and use the power spectrum 
and the density of states values at this energy: S x (u>) = S x (E n ,uj) and p = p(E n ). For 
the calculation of the absorption coefficient the global dependence on the energy must be 
restituted. 

The power spectrum S x (u) has the symmetry property S x {—uj) = S x (uj). The depen- 
dence S x (uj) in the domain uj > in typical strongly chaotic systems, such as nonlinear 
oscillators [20] and billiards [21, 22], has the form of an asymmetric peak. We shall define 
the peak value of the Rabi frequency simply as Q, the frequency of the maximum as Co and 
the characteristic width of the peak as T. Typically the ratio uj/T is about few units. 

Right after the switching on the perturbation all amplitudes (but that of the initially 
populated state) grows in absolute value linearly in time. At this ballistic stage the energy 
dispersion grows quadratically in time: 

AE 2 ) w K^W^Tpt 2 , (10) 

where K\ is a numerical constant. This stage is limited by the depletion of the initial 
population and lasts until the depletion time ~ (ftFp) -1 / 2 . At this time considerably 
populated levels are spread over the energy range AE ~ huj that contains many levels (since 
p oc TT d with d > 2). We shall expect that at the next stage the ensemble averaged density 
of probability is a smooth function with a characteristic scale AE » hu. 



It is convenient to write indices in Eq. (9) as arguments of functions. We shall use the 
frequency distance from the initial level as a basic independent variable u>, cik will be denoted 
as a(e), where e = (E k — E n )/h. Dummy variables rj and rj will have the same meaning. By 
formal integration of Eq. (9) and subsequent recurrent substitution we obtain the equation 
for the rate of change of the local probability density w (e) = \a {s)\ 2 : 



dw (e) _ ^ n ^ ^ ei{e _ v)t ^ fl t ) x (n) 



dt 

t 



x J dt'tt (e, n') e~ iie ~ n,)t cos u t a* {rj , t) + c.c. 



Summation in this formula goes over discrete values of rj and rj', and this equation is still 
exact. 

Now we assume, that the amplitudes a (rj, t) are random processes, that are not correlated 
for different states: (a(rj,t) a* (?/,£')) ^ <W- Then for the averaged probability density we 
can retain in the Eq. (11) only diagonal terms: 

d(w(e)) 



dt 



^2(n 2 (e,r])\ cosu tx (12) 



x / dt' e ^-^*-*') coscj t' (a (rj, t) a* (rj, t')) + 



c.c. 



With the assumption that the averaged (w(e,t)) is a smooth function of e, and slowly 
varying function of t, we may cast the product of amplitudes in the form 

(a ( V , t) a* (rj, t')} + c.c. = 2 (w (rj, t)) B (t - t') , (13) 

where B (r) is the normalized (B (0) = 1) autocorrelation function of the probability ampli- 
tudes. 

By replacing the averaged square Rabi frequency by its value from Eq. (3) (that depends 
only on the difference e — rj), substituting the summation over the states by the integration 
weighted with the density of states, and averaging over the time intervals that are much 
larger than the field period, we come to the equation 



dw (e) 



= J drfVt 2 (e — rj) p (rj) J dr cos (e — rj) t cos lu tB (t)w (rj) . (14) 



dt 

o 



From now we drop the angular brackets and deal only with ensemble averaged quantities. 
If the rate of variations of w (rj, t) is small in comparison with the decay of correlations of 



amplitudes given by B{r), then we can treat Eq. (13) as summation over the probability 
flow that comes from the different parts of the frequency range with the constant rate, 

oo 

W (77 — > e) = Vt 2 (e — rj) J dr cos (e — 77) r cosluotB (r) . (15) 



This approximate expression to some extent replaces the Fermi golden rule for strong per- 
turbations. 

To construct the kinetic equation one must take into account both incoming and outgoing 
flows of probability. By taking into account the total probability flow, expanding w (e) in 
the Taylor series, we obtain a diffusion equation with the probability diffusion coefficient in 
the energy scale 

00 

-D ~ J drjh 3 ri 2 tt 2 (ri)pJ(ri,oj ), (16) 

—00 

where 

00 

J(r],co )= J dr cos rjr cos uj tB (t) . (17) 


Now the problem reduces to the calculation of the integral J(r),cuo). For large enough 
times the average probability density, that is governed by the diffusion equation, varies 
slowly, and we can treat the system Eq. (9) as a set of equations in which all a m (t) are 
non-correlated random processes with the same statistical properties. Then by averaging 
the equation for the squared time derivative of an amplitude, we obtain the expression for 
the mean squared frequency of these processes, that also gives the estimate for the square 
of the correlation decay rate 7: 

(u> 2 ) = ya i (Ti)hpdT, = >r 2 . (is) 

From Eq. (18) for the decay correlation rate we have the estimate 



7^X2^0^, (19) 

where K 2 is a numerical constant In the domain of strong field the autocorrelation function 
is the fastest component under the integration sign in Eq. (17). Then we have J ~ 7 _1 ~ 
. By the substitution of this expression into Eq. (16) we obtain the estimate 
of the energy diffusion coefficient: 

D w Kz% 2 Cj 2 VL^rhp, (20) 

where K3 is a numerical constant; this expression agrees with the dependence that was 
derived in Sec. 1. One must recall that this expression is valid only for the nearly resonant 



perturbation frequency ( \u — cu\ < T ). We do not enter here into the studies of the 
dependence of D on the perturbation frequency ujq in the wider domain, postponing it for 
the further studies. 

III. NUMERICAL EXPERIMENT 

To check the analytical results of the preceding section, the system of equations (9) 
has been integrated numerically. The number of equations varied from iV = 300 to N = 
1200 with the purpose to suppress the influence of the border. The envelope of the Rabi 
frequencies has been taken in the double lorentzian form 

r r 2 r 2 l 

Qm,n = ^ - ~2 . -po 7 , ~\2 , -P2 ' 

_{Lo m -uj n + u) +T 2 (u m -u n + u) +T\ 

All calculations have been carried for the "resonant" perturbation frequency uo = lj and for 
the peak width T = 0.3a). 
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FIG. 1: The dependence of the logarithm of the probability density on the dimensionless frequency 
e' = ejuo for the time values t = 5a) _1 (a), t = 10u) -1 (b), and t = 15cD _1 (c). The grassy lines - 
values of \nw (e'), averaged over 10 different sets of matrix elements, solid lines - fitted parabolas. 
To avoid the overlap of graphs, the plots for cases (b) and (c) are shifted upwards by 3 and 6 units 
respectively. 



Fig. 1 shows the distribution of probability as a function of dimensionless frequency 
e' = s/uj for different moments of time. It is clearly seen that even for relatively small 
time t = 5a; -1 = 2.2t rf the distribution has very accurate Gaussian form, with deviations 
noticeable only for \e'\ > 7.5. Thus we have quantitative support for our conclusion about 
the diffusive character of the energy evolution. 
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FIG. 2: Dependence of the autocorrelation function of the probability amplitudes on the dimen- 
sionless time shift r' = tuj for three different sets of parameters with uhp = 30 (black squares), 
ujhp = 60 (open circles) and tohp = 120 (black triangles) and the same value of Qy/Thp = 0.618a). 
The statistical errors are about the size of the data symbols. 

Fig 2 depicts the form of the normalized autocorrelation function of the probability 
amplitudes. Values B (r) have been calculated numerically for three sets of parameters with 
different values of p but the same value of the product Q^/KTp. It is clearly seen that B (r) 
for these sets are nearly identical, as expected. The decay rate 7 taken from the equation 
B (7^) = exp (—1) is 7 = 0.77 VtJhTp , that supports the estimate Eq. (19). 
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FIG. 3: The dependence of the ratio R = D/2Wf of the energy diffusion constant D to the doubled 
Fermi transition rate Wf = (7r/2)Q 2 hp on the logarithm of the ratio of the Rabi frequency to its 
border value L = fi/fij,, = (2tD /nhp) 1 / 2 . The dashed line represents the curve R = Aexp (-L) 
that corresponds to the theoretical dependence Eq. (20); it is fitted to the last three points. 

Figure 3 represents the dependence of the ratio R = D^2Wp of the energy diffusion 
constant D to the doubled Fermi transition rate W F = (n/2) il 2 hp on the logarithm of the 

— 1 /2 

ratio of the Rabi frequency to its border value L = Q/Qj,, = {2u)/iiTip) . It is seen that 
for the weak field this ratio comes close to the asymptotic limit (unity), decreases in the 
vicinity of the border and decays as F -1 for sufficiently strong fields. The agreement with 
the theoretical estimates is quite convincing. 

IV. CONCLUSION 

From the comparison of the numerical data with the theoretical estimates we can con- 
clude that the approach of the Sec. 2 gives reasonably accurate description of the energy 
evolution process in strong fields, in spite of numerous simplifying approximations used in 
the calculation. 

To improve the accuracy and, firstly, to derive from the first principles the equation for 



the correlation function B(r), the model of the random process must be improved. We used 
the model of the stationary process, whereas from Eq. (9) one can conclude that the model 
of periodic random process would be more appropriate. 

The main conclusion from the results of our calculation is a qualitative one: by substi- 
tuting Eq. (3) in Eq. (16), we obtain D xs K 3 huj 2 F [S x (u> ) T] 1 ^ 2 . This quantity in the 
classical limit h — > vanishes along with the energy absorption rate Q (see Eq. (5)). That 
means the violation of the quantum-classical correspondence for the absorption - and, more 
generally, for the linear susceptibility of chaotic systems to harmonic external fields. 
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